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Abstract. Few-nucleon correlations in nuclear matter at finite densities and 
'sj" i temperatures are explored. Using the Dyson equation approach leads to effec- 

I tive few-body equations that include self energy corrections and Pauli block- 

■ ing factors in a systematic way. Examples given are the nucleon deuteron 

in-medium reaction rates, few-body bound states including the a-particle, 
^ and a-particle condensation. 

> . 

1— I ' 1 Introduction 

O 

pH I Strongly correlated many-particle systems provide an exciting field for applica- 

tions of few-body methods. Examples are nuclear matter, quark matter (quark 
gluon plasma), ionic plasmas, among others. In an equilibrium situation at fi- 
nite temperatures and densities one usually introduces a mean field to describe 
' the gross features of many-particle systems. However, many exciting phenom- 

ena, such as the formation of clusters below a certain density (Mott density) and 
the appearance of superconductivity below a critical temperature, cannot be de- 
scribed in a framework of noninteracting quasiparticles (ideal system). Within 
a quantum statistical approach it is possible to go beyond such a picture and 
consider the residual interactions between the quasiparticles. A well known equa- 
tion, e.g., to study two-body correlations is known as Bethe-Goldstone equation 
and/or Feynman-Galitskii equation, depending on details Q. To go beyond two- 
body correlations is desirable for many reasons. Among them are: i) For a mi- 
croscopic description of the heavy ion collision, three-body reaction rates are 
an important input into the collision integral. Their medium dependence has 
hardly been studied |^, ^, 0, ^. ii) Bound states are effected by the density and 
temperature of the medium and the binding energy may become zero leading to 
the Mott effect |, |, |, |, 0. Hi) Since the a-particle is the strongest bound 
nucleus it should be relevant for the equation of state of nuclear matter. It might 
induce a-condensation and quartetting which would be a different state of mat- 
ter besides superfluidity induced by pairing or pair condensation [11]. iv) The 
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three-quark system is particularly interesting because it is a reasonable model of 
baryons and useful to study the phase transition between the quark phase and 
the hadronic phase of nuclear matter |12, 13|. v) The three-body input into the 
two-body t-matrix implied by the equations hierarchy has not been studied and 
might, e.g., affect the critical temperature. 

The first three issues will be addressed in the following, iv) is given to some 
detail in Ref. [^] of this issue, and v) is left for future investigations and is 
relevant, e.g., for the question of a possible color superconducting phase. 

2 Theoretical tools 

We utilize the Dyson equations to tackle the many-particle problem. An review is 
given in Ref. [14|. This enables us to decouple the hierarchy of equations finally 



deriving effective in-medium few-body equations. The many-particle Hamilto- 
nian is given by 

H = Y^ Ho{l)a\ai + ^ Yl ^2(12, l'2') 44a2.ar (1) 

1 121'2' 

where Hq is the single particle energy and V2 a generic two-body potential. The 
chronological equal time cluster Green function for a given number of particles 
is defined by 

or/ = -i{TMT)Al{T')), (...) = Tr{po...} (2) 

where Aa = 010203 . . . o„, e.g., are Heisenberg operators and pQ the equilibrium 
density operator. The respective Dyson equation for the cluster is 



i^Or/ = S{t-t'){[A^,aIU)+J2 [dfMl-^g;/. (3) 



d_ 



The mass matrix that appears in ^ is given by 

Ml/ = 5{r - T')Mlp^, + (4) 

cluster mean field retardation 
{MM)lp^, = {[[A^,H]{tIAI{t)]) (5) 

{MMy[-/p = Y.^Tr[A^,H]rAA\,HU),,,,^^,^,,. (6) 



More details are given in [14|. The equation (^) is expressed in momentum space 
and the time component as Matsubara frequency r — > that is analytically 
continued into the complex plane, zx z, for a textbook treatment see 
To arrive at suitable calculable expressions the following approximations are 
utilized: i) Only the cluster mean field contribution to the kernel (Q) is used; ii) 
The density operator p is evaluated for an uncorrelated medium. This way the 
equations hierarchy is decoupled and effective few-body equations that describe 
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few-body correlations including medium effects have been derived. For = ai, 
the one-particle Green functions is 

G{z) = {z - e^)-' (7) 

where the quasi-particle self energy is 

The last equation introduces the effective mass that is a valid concept for the 
rather low densities considered here and ^'^^ ^ 1^ — ^^^(0). The Fermi function 
fi = f{^i) for the i-th particle is given by 

The resolvent Go for n noninteracting quasiparticles is 

n 

Go{z) = {z-Ho)-^N = Ro{z)N, Ho = J2^i (10) 

i=l 

where Go, Hq, and are formally matrices in n particle space. The Pauli- 
blocking for re-particles is 

N = flf2...fn±flf2...fn, / = !-/ (H) 

where the upper sign is for Fermi-type and the lower for Bose type clusters. Note 
that NRq = RqN. Straight forward evaluation of (^) using the Wick theorem 
leads also to the full resolvent G{z). For each number of particles in the cluster 
the resolvents have the same formal structure and can be written in a convenient 
way close to the one for the isolated system, viz. 

G{z) = (z-Ho- Vy'N ^ R{z)N, ^ ^ ^2°^2"- (12) 

pairs a 

Note that ^ V and R{z)N ^ NR{z). To be specific, for an interaction in 
pair a = (12) the effective potential reads 

(123 . . . liv(^'V2('')|l'2'3' . . .) = (/1/2 - /i/2)V2(12, 1'2')<533' • • • (13) 



For further use in the Alt Grassberger Sandhas (AGS) equations [15] we give 
also the channel resolvent 

G^z) = {z-Ho- N^V^^'y'N = R^{z)N. (14) 

For the two-body case as well as for a two-body subsystem embedded in the 
re-body cluster the standard definition of the t matrix leads to the Feynman- 
Galitskii equation for finite temperature and densities [^, 

T^{z) = V2 + V^N!^Ro{z)T^{z). (15) 



4 Three- and Four-body correlations 

Introducing the AGS transition operator via 

R{z) = 5apRa{z) + Ra{z)Ua^{z)Rp{z) (16) 
the effective inhomogeneous in-medium AGS equation reads 

Uap{z) = (1 - 5ap)R^Hz) + N]T^iz)Ro{z)U^piz). (17) 

The homogeneous in-medium AGS equation uses the form factors defined by 

7 

to calculate the bound state V'-Bs- Because of the non-symmetric form of the 
potential the equation for the form factors and the dual are different 

\Fa) = Y.~^a0N^T^{B^)Ro{B^)\Fp) (19) 
\Fa) = Y.~^apT^{B^)N^R^{B^)\Fp). (20) 

Finally, the four-body bound state is described by 

= 2^^- U^.^{Bi) Ro{Bi)N] T^{Bi) Ro{B^)\J^), a C a,j C r. (21) 
""^ 3-body 2-body 

where u, r denote the four-body partitions. The two-body input is given in ([l5| ) 
and the three-body input by (p!7|), both medium dependent. 

3 Results 

3.1 Reaction rates 

An experiment to explore the equation of state of nuclear matter is heavy ion 
collisions at various energies. Here we focus on intermediate to low scattering 
energies and compare results to a recent experiment ^■^^Xe-|-^^^Sn at 50 MeV/A 
by the INDRA collaboration |jl^. A microscopic approach to tackle the heavy ion 
collision is given by the Boltzmann equation for different particle distributions 



F = ifp, fn, Id, ft, fh) up to /i and t 

^^F{t) + {e, Fit)} = ICnFit)] {Fit) ± 1) - )Cr[F{t)] Fit) (22) 

where e is the mean field energy and IC denotes the respective collision integrals 
that include, e.g., the one for deuteron loss Nd NNN, 

ICT{P,t) = Jd^kJ d^hd^k2d^h \{kik2ks\Uo\kP)\l^^p^j^ 

medium dependent! 

xfNikl,t)fNik2,t)fNik3,t)fNik,t) + ... (23) 
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where / = 1 — /. A solution is given via a Boltzmann Uehling Uhlenbeck (BUU) 
simulation |l^ . As indicated in ( p3|) the reaction rate is in principle medium 
dependent. However, previously this medium dependence has been neglected. 
Within linear response theory for infinite nuclear matter the use of in-medium 
rates leads to faster time scales for the deuteron life time and the chemical 
relaxation time as has been shown in detail in Refs. ^. 
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Figure 1. BUU simulation of the deuteron 
formation during the central collision of 
^2^Xe-K"^Sn at 50 MeV/A. 
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Figure 2. Ratio of proton to deuteron num- 
bers as a function of cm. energy. The exper- 
imental data are from the INDRA collabora- 
tion. 



Now we use the in-medium AGS equations (|r^) that reproduce the exper- 
imental data in the limit of an isolated three-body system. For details on the 
specific interaction model see Ref. Q. We investigate the influence of medium 
dependent rates in the BUU simulation of the heavy ion collision as compared 
to use of isolated (i.e. experimental) rates. Figure |^ shows that the net effect 
(gain-minus- loss, eq.(|2^)) of deuteron production becomes larger for the use of 
in-medium rates (solid line) compared to using the isolated rates (dashed line). 
The change is significant, however, a comparison with experimental data is dif- 
ficult since deuterons may also be evaporating from larger clusters that has not 
been taken into account in the present calculation so far. The ratio of protons to 
deuterons may be better suited for a comparison to experiments that is shown 
in Figure |. The use of in-medium rates (solid line) lead to a shape closer to the 
experimental data (dots) than the use of isolated rates (dashed line). 



3.2 Bound states, Mott effect 

In these calculation, besides the change of rates, also the Mott effect has been 
taken into account. Figure ^ shows the dependence of the binding energy for 
different clusters at a given temperature of T = 10 MeV and at rest in the 



medium. At first sight an Efimov effect [17| might be expected in the vicinity 
of the Mott transition of the deuteron. However, two main reasons prevent the 
Efimov effect to appear in a simple way: i) The deuteron binding energy in 
the medium depends parametric (through the blocking factors) on the deuteron 
momentum. Since the deuteron-like subsystem is not at rest, in other words the 
effective strength of the potential that enters into the three-body problem varies 
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Figure 3. Difference between tfie pole en- 
ergy of the bound state and the continuum, 

B{n, T) = Epole — -Bcont- 
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Figure 4. Critical temperatures of condensa- 
tion/pairing leading to superfluid nuclear mat- 
ter. For an explanation see text. 



with momenta, a possible Efimov effect is washed out. ii) The excited states that 
should appear from the continuum (Efimov states) are as well blocked by the 
medium. This blocking may not be so strong as the ground states, because the 
momentum distribution is peaked at higher momenta. As a consequence only a 
careful quantitative analysis might answer the question of Efimov states. On the 
other hand nuclear matter might not be the best system to eventually observe 
such an effect. 

In Figure ^ part of the phase diagram of nuclear matter is shown. The lines 
indicate phase transitions. The critical temperatures of condensation/pairing 
(dashed line, ||ll|) leading to superfluid nuclear matter are shown. The possible 



area of a condensation (solid line) as suggested by |11] is also given. The latter 
is based on a variational calculation using the 2-1-2 component of the a parti- 
cle to evaluate the condition for the onset of superfluidity for the four-particle 
system B(Tc,fi,P = 0) = 4/x. The critical temperature found by solving the 
homogeneous AGS equation for /x < confirms the onset of a condensation even 
at higher values (dotted line). For /i > the condition B = Afi for the phase 
transition can also be fulfilled. However, the homogeneous AGS equation cannot 
be used to investigate the steep fall-off predicted in Ref. because of contin- 
uum poles that are not compensated by the blocking factors of the potential as 
is the case for the two-body problem. Whether the steep fall-off is of physical 
origin or due to the use of a homogeneous equation also for /x > needs further 
investigation. 
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